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We have measured the thermal conductivity of the iron pnictide superconductor LaFePO down to 
temperatures as low as T=60mK and in magnetic fields up to 5 T. The data shows a large residual 
contribution that is linear in temperature, consistent with the presence of low energy electronic 
quasiparticles. We interpret the magnitude of the linear term, as well as the field and temperature 
dependence of thermal transport in several pairing scenarios. The presence of an unusual supralinear 
temperature dependence of the electronic thermal conductivity in zero magnetic field, and a high 
scattering rate with minimal T c suppression argues for a sign-changing nodal s± state. 



INTRODUCTION 

One of the most important aspects of formulating a de- 
tailed and complete understanding of superconductivity 
in the iron-based superconductors is knowing the sym- 
metry of the superconducting gap, because of its impli- 
cations for the pairing mechanism. In the presence of 
a weak electron-phonon interaction [T], focus has cen- 
tered on a pairing mechanism associated with magnetic 
fluctuations [2]. The combination of a magnetic cou- 
pling mechanism and the multiband nature of the Fermi 
surface gives rise to the possibility of different sheets of 
the Fermi surface with superconducting gaps of different 
phase and symmetry. The permutations of these possi- 
bilities leads to a very rich phenomenology, and makes 
the interpretation of experimental results far more diffi- 
cult than in a single band situation, as for example in 
the cuprates. Moreover, theoretical calculations [5] have 
shown that many of these potential gap arrangements lie 
very close in energy and therefore each requires careful 
consideration. 

Experimental studies across the various members of 
the pnictide family of superconductors have revealed a 
complex picture [I]. NMR Knight shift measurements 
establish that the superconductivity in the iron pnic- 
tides is singlet in nature [SH5] , and penetration depth 
measurements [9], thermal conductivity [10] , Andreev- 
reflection spectroscopy [11] and Angle Resolved Pho- 
toemission (ARPES) studies [12], have revealed fully 
gapped superconductivity in many compounds. How- 
ever low temperature penetration depth measurements 
in BaFe 2 (Asi_ x P x ) 2 [13] and LaFePO [Hj, and Raman 
scattering in LaFePO [15] . observe low energy excitations 
suggestive of a line node. Attempts to resolve these ob- 
servations with theory have highlighted the possible role 



of intraorbital interactions [TH [17] and the precise details 
of Fermi surface nesting |18) in determining different gap 
topologies. It has even been shown that in the case of 
an isotropic sign-changing s-wave (s±) state, under cer- 
tain conditions disorder can create subgap states which 
can give rise to a low energy density of states [19] mim- 
icking the presence of nodes. The closeness in energy of 
these various gap structures suggests that different ma- 
terials may support different gap topologies, which en- 
sures a rich phenomenology across the pnictide materials. 
Clearly, it is extremely important to conduct high-quality 
experiments in many compounds in order to establish 
whether the low-energy excitations are intrinsic (arising 
from nodes) or extrinsic (arising from disorder) so that 
scenarios for superconductivity in this multi-band system 
may be tested. 



To shed light on this issue, we use thermal conductiv- 
ity measurements as a bulk probe of the superconducting 
state. By carrying out our studies at temperatures as low 
as 40 mK and in the presence of an applied magnetic field, 
we are able to accurately separate out contributions from 
electrons (n e ) and phonons (K p h)- Our study unambigu- 
ously reveals the presence of a superconducting gap with 
line nodes, by establishing a non-zero linear component 
of thermal conductivity as the temperature approaches 
absolute zero, which arises from low-energy quasiparti- 
cle excitations. A quantitative analysis of the magnitude 
of the electronic thermal conductivity and its evolution 
with temperature and applied magnetic field is used to 
attempt to discriminate between d-wave and nodal s± 
superconductors. 
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EXPERIMENTAL 

Our study focuses on LaFePO, an iron-phosphide 
superconductor with T c ~ 7.5 K, and a reported 
H C 2 ~ 900mT [20], which can be made with very 
high purity. LaFePO is nonmagnetic, isostructural with 
LaFeAsO and has a Fermi surface consisting of almost 
nested electron and hole pockets [5T]. Our samples were 
small, single crystal platelets of dimensions ~ 200 /im x 
50 /im x 20 /im grown via a Sn flux method [22] . Electri- 
cal resistivity measurements confirmed the high quality 
of our samples, revealing a residual resistivity ratio of 
P300k/pok = 65 with p = 2.4 /iJ7cm, measured by sup- 
pressing superconductivity with H — 2 T applied along 
the c-axis. 

Contacts with very low thermal resistance are essential 
for measurements of thermal conductivity at T < IK, 
as thermal decoupling between the electron and phonon 
degrees of freedom is known to potentially lead to an un- 
derestimate of K e [23] . Highly conductive contacts were 
prepared by Ar etching the crystal surface and evaporat- 
ing Pt pads, which were bonded to gold wires with silver 
epoxy. At low temperatures the electrical resistance of 
these is less than 1 m£l. 

The thermal conductivity was measured using a single 
heater-two thermometer method. The heat current was 
supplied along the ab planar direction direction, and the 
magnetic field applied perpendicular to this. The mea- 
surements were made in a dilution refrigerator by varying 
the temperature from 0.04 K to > 0.7 K at fixed mag- 
netic field. To ensure homogeneous flux penetration, the 
samples were field-cooled by cycling to T > T c before 
changing the field. The error in the absolute value of 
the conductivity is estimated to be approximately 10%, 
which is set by uncertainties in the geometric factor of the 
sample. The relative error between temperature sweeps 
at different fields however is lower, on the order of 3%. 



RESULTS AND DISCUSSION 

Fig. [T] shows the raw measurements of the in-plane 
thermal conductivity plotted as k/T versus average sam- 
ple temperature (T avg ) for temperatures from 60 mK to 
1.0 K, for a variety of magnetic fields applied parallel to 
the c-axis and perpendicular to the heat current direc- 
tion. In a simple non-magnetic metal, the temperature 
dependence of the thermal conductivity is expected to be 
governed by two terms - an electronic linear term n e cx T 
and a cubic phonon term n p h cx T 3 . In the superconduct- 
ing state, K e is dramatically altered as it is determined 
solely by thermally excited quasiparticles resulting from 
broken Cooper pairs. Below T c , the low-energy landscape 
of excitations depends sensitively on the topology of the 
superconducting gap. 
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FIG. 1: Temperature dependence of thermal conductivity di- 
vided by temperature at magnetic fields applied perpendicular 
to the c-axis. Inset: Thermal conductivity divided by tem- 
perature and charge conductivity in thermal units (Lo/p) at 
H=2 T as a function of temperature. The agreement between 
the two independent measurements shows that the Wiede- 
mann Franz law is satisfied. 



Beginning with the zero-field data, we note that the de- 
pendence of k/T on temperature is much stronger than 
any of the in-field data, which suggests an extra tem- 
perature dependent contribution to thermal conductivity 
that is rapidly suppressed with field. This qualitative fea- 
ture has been seen in the superconducting states of both 
cuprate [2J and filled skutterudite [25] superconductors. 
Since the phonon term is not expected to be particularly 
sensitive to field, we interpret this extra zero-field con- 
ductivity as arising from thermally excited quasiparticles, 
in addition to the expected linear term. 

We now set out to separate the phonon and electron 
contributions to k. When a field of 20 mT or greater is 
applied, the temperature dependence of k/T is greatly 
reduced, and remains unchanged for higher fields. We 
thus assume that for H > 20 mT, phonons are the only 
other contribution to the conductivity in addition to the 
linear electronic term. The phonon term is expected to 
be small at low temperatures, as demonstrated by the 
relatively small slopes of the data for H > 20mT in Fig. 
[T] We then take as our phonon contribution the total 
conductivity measured at H — 20 mT less the residual 
linear term kq/T at that field (we define kq/T as the 
limit of k/T as T -> 0). 



K p h k(20 mT) 
~T~ ~ f 



k (20 mT) 
T 



(1) 
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This phonon conductivity is plotted in Fig. [2] We 
assess the magnitude of K p h by fitting to a cubic temper- 
ature dependence consistent with scattering of phonons 
at the boundary of the sample, giving K p h = 1.2 T 3 
mW/K 4 cm. Theoretically, the coefficient is given by ki- 
netic theory to be equal to the product of the phonon 
heat capacity, /3, the sound velocity, v s and a mean-free 
path equivalent to the geometric average of size of the 
sample cross-section (a x b), I 



Kph 



ph 



-f3T 3 v s l ph 



(2) 



Measurements of specific heat [21] give a phonon con- 
tribution of /3 = 0.16 mJ/K 4 mol and we can estimate an 
average sound velocity of v s — 6263 m/s from the De- 
bye temperature, Qjj = 371 K. Combining this with our 
sample dimensions l p h — 2.7 x 1CP 5 m, gives a phonon 
conductivity of n p h = 1 T 3 mW/K cm which is consis- 
tent with our measured value. 

We can now use this phonon conductivity to extract 
the entire electronic conductivity in zero field by taking 



K e (H = 0) k(H = 0) 



T 



T 



(3) 



This yields the n e /T curve labeled 'Electronic' shown 
in Fig. [2] which evolves linearly with temperature be- 
tween 60 mK and 600 mK, making the total electronic 
conductivity the sum of a linear and quadratic term. 
When the the data is extrapolated to zero temperature, 
the residual linear term is hq/T — 3.0 ± 0.3 mW/K 2 cm. 
Extrapolating the raw total conductivity fitted in the 
same way over the same temperature range leads to a 
linear term kq/T = 2.9 ± 0.3 mW/K 2 cm, indicating the 
conclusion is not significantly dependent on any analy- 
sis. The magnitude of the normal state conductivity is 
taken as the linear term at 700 mT, k (0.7 T)/T = 9.6 
± 0.3 mW/K 2 cm. Consequently, the magnitude of the 
zero-field residual linear term is a substantial (0.31) frac- 
tion of the normal state conductivity. We note that our 
current study extends thermal conductivity data to al- 
most an order of magnitude lower in temperature than 
previous work [20., indicated by the filled circles in Fig. 
[2] Although there has been no attempt to subtract any 
phonon conductivity, the data is qualitatively and quan- 
titatively consistent with that reported here. This allows 
us to make firmer conclusions about the magnitude of 
k /T, to capture any temperature dependence of K e /T, 
and to rule out the opening of a low energy gap between 
60 and 600 mK. 

Having extracted n e /T(H = 0), we may now com- 
pare our data to quantitative predictions for the size and 
temperature dependence of the electronic conductivity 
expected to arise for various gap symmetries [2TJ [28] . 

Fully gapped s± symmetry: Conventional thinking re- 
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FIG. 2: Thermal conductivity divided by temperature at zero 
magnetic field. We show the total thermal conductivity, as 
well as the estimated electronic and phonon contribution and 
the calculated electronic contribution in the d-wave Born and 
unitary scattering limits (described in the text) . Closed circles 
are data from [20l. 



garding a fully-gapped superconducting state leads to 
the conclusion that the residual electronic conductivity 
is zero because the quasiparticle excitations are exponen- 
tially suppressed in the low temperature limit. Interest- 
ingly, in the sign changing s± state, where isotropic gaps 
of different signs exist on different Fermi surface sheets, 
interband scattering can lead to a small but finite den- 
sity of states even at zero energy |27j . The magnitude 
of the thermal conductivity associated with this band re- 
mains below 1 % of the normal state conductivity though, 
even for sizable scattering rates. Consequently the data 
reported here for LaFePO is quantitatively inconsistent 
with a fully-gapped s± symmetry. 

Nodal s± symmetry: Depending on the exact details 
of the magnetic fluctuations that provide the electronic 
coupling, the s± symmetry can develop minima on one 
of the Fermi surfaces. These minima can either be very 
deep, retaining the same sign of the order parameter, or 
may actually change signs. Many transport models treat 
the order parameter in an effective two band model, with 
an isotropic gap on the sheet centered at the T point, 
with the sheet about the M point playing host to the 
anisotropic Fermi surface [3]. 

As firmly established in a d-wave superconductor, a su- 
perconducting gap with nodes in the presence of impurity 
pair-breaking leads to an impurity band of quasiparticle 
excitations that give rise to a residual linear term in the 
electronic thermal conductivity. The essential difference 
between the conductivity in the nodal s± case and the 
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d-wave case is that, in the latter, the conductivity is uni- 
versal in the sense that the magnitude does not depend 
on the scattering rate. This is not the case for the s± 
symmetry as shown in calculations by Mishra et al. [57] 
where the magnitude of the conductivity is strongly de- 
pendent on the normal state scattering rate. 

The normal state scattering rate for the sample re- 
ported here can be estimated from the normal state elec- 
tronic conductivity, K e (iV), using kinetic theory which 
gives K e (N)=jv F r/3. The average magnitude of the 
electronic specific heat 7 measured by two groups is 
C e /T = 7 ~ 11.5 mJ/K 2 mol [26l ES] and the Fermi 
velocity is vp — 1.5 x 10 5 m/s from ARPES measure- 
ments [30 . This leads to a normal state scattering rate 
r = l/2r = 1.1 x 10 12 s _1 . In reduced units with respect 
to T c , this is fiT/k B T c = 1.1. 

An alternate method is to estimate T from de Haas- van 
Alphen measurements of the mean free path £q. Using £q 
« 1000 A [22 and v F « 1.5 x 10 5 m/s [30. yields r = 
vf/£o = 1 x 10 -13 s. This would suggest our samples of 
LaFePO are in the regime TiT /ksT c ~ 0.8, consistent with 
the previous method. Similar values have been observed 
in other pnictide compounds, for instance hT/ksTc ~ 
1-2 has been inferred from transport and penetration 
depth measurements in Ba(Fei_ K Coa;)2As2 |10| . LiFeAs 
[3T] and Bai_ ;E K :i; Fe2As2 [32]. 

This high value of scattering rate might be expected 
to result in significant T c suppression, however this does 
not appear to be the case in LaFePO. Indeed T c has been 
shown to be relatively insensitive to disorder, increasing 
Po by a factor of 5 has no noticeable effect on T c [33]. 
These observations are consistent with theoretical models 
describing how nodes nodes of an extended-s state may 
be lifted by disorder, leading to a disorder induced T c 
suppression much more gradual than that expected from 
Abrikosov-Gor'kov theory |34) . 

In both the 'deep minima' and sign changing node sce- 
nario, a scattering rate of this order of magnitude can 
yield a value of kq /T that is a sizable fraction of the nor- 
mal state value [27], although other parameters such as 
the relative size and number of the gap maxima are re- 
quired for a detailed quantitative comparison. At finite 
temperatures, the extra T 2 contribution to K e is more 
consistent with the sign changing nodes scenario. For a 
scattering rate of hT/kBT c ^0.1, assuming similar sized 
gap maxima on both sheets, n e /T is predicted to behave 
approximately linearly up to T=0.2 T c , as in our data. 
For a reasonable range of parameters in the 'deep min- 
ima' scenario, a constant value of n e /T is expected [27], 
which is not reflected in our data. A detailed calcula- 
tion using measured parameters for LaFePO and a full 
account of the details of inter versus intraband scattering 
would be useful in distinguishing between these scenarios. 

d-wave symmetry: The thermal conductivity of a d- 
wave superconductor has been studied in great detail 
both experimentally and theoretically. Theory predicted 



a universal linear electronic conductivity with a second 
order term that is cubic in temperature and depends 
strongly on scattering rate and strength [35] . 
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where 7 is the impurity bandwidth and a is a constant 
that depends on the scattering strength. Both these at- 
tributes have been observed in high-temperature cuprate 
superconductors [HJ [35] . In this case, we can estimate 
the value of the universal linear term for LaFePO using 



T 



4 nr 1 \ n n 



(5) 



Taking the weak coupling <i-wave gap expression, Ao = 
2.l4ksT c , and that fi — 2 represents the slope of the 
gap at the node for a pure d-wave function, we compute 
the universal linear conductivity, n /T = 2.7 mW/K 2 cm. 
This is in excellent agreement with the measured value 
shown in Fig. [2] Turning to the second term in Eq. [4] 
the main observation is that the temperature dependence 
is linear, i.e. the low temperature correction to the linear 
term is quadratic not cubic. To better evaluate the level 
of consistency, we can estimate the magnitude of the con- 
ductivity in this model using the normal state scattering 
rate in both the limit of strong (unitary) and weak (Born) 
scattering. In the unitary limit, the impurity bandwidth 
is given by 7 = 0.63\/A^r = 8.5 K and a = 1/2 [55] . 
This gives the curve shown in Fig. [2] for unitary scatter- 
ing. On the same plot, we show the expectation for the 
weak scattering limit computed for a = (ttv2Tq)/2 and 
7 = 4A exp(-7rA /2r) = 0.1 K. 

Clearly the magnitude of the measured conductivity is 
bounded between these two limits indicating the possi- 
bility that the scattering phase shift lies between and 
7r/2, and consequently the range of applicability of the 
low temperature correction as defined by 7 will lie be- 
tween 0.1 - 8.5 K. The possibility of anisotropic scat- 
tering should also be considered as a means to recon- 
cile experiment and theory in this case. Consequently, 
a comparison with the d-wave symmetry model provides 
quantitative agreement for the linear term on one hand, 
but a qualitative (not cubic) and quantitative discrep- 
ancy with the temperature dependence of the electronic 
conductivity above the linear term. 

An important test of the various gap scenarios dis- 
cussed above is investigating the sensitivity of n e /T to 
scattering rate, which is expected to be negligible in the 
d-wave universal limit and be substantial in the s± sce- 
narios. In the previously reported LaFePO study, we note 
that the normal state conductivity is 6.0 mW/K 2 cm [20] . 
compared to 9.6 mW/K 2 cm in the present study. This 
suggests that the scattering rate for the previous study is 
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FIG. 3: Main panel (colour online): The linear, electronic con- 
tribution to thermal conductivity normalized to the value in 
the field-induced normal state. The data is plotted as a func- 
tion of magnetic field in units of £f C 2, with the blue squares 
showing the extrapolated T=0 value, and the open circles 
showing the T=0.46 K value. The closed circles are taken 
from ref [20]. The inset shows a zoom of the low field region. 



1.6 times higher than for this work, however the extrap- 
olated linear term in both cases are ~ 3.0 mW/K 2 cm. 
Within a d-wave scenario, the prospect of a universal 
conductivity is not ruled out assuming a linear extrap- 
olation and an error in the absolute magnitude of each 
measurement of 10% due to measurements of the geo- 
metric factor. A controlled study looking at different lev- 
els of disorder in this material would prove very useful in 
confirming universality and thus distinguishing between 
the d and s± scenarios. 

The application of a magnetic field can be used to tune 
a superconductor between a fully superconducting and a 
fully normal state, which can provide valuable informa- 
tion about the gap topology via the contribution of field- 
induced quasiparticles to thermal transport. The inset 
to Fig. [T] shows the normal state thermal conductivity 
divided by temperature, accessed by applying H — 2T 
/ /c. Also plotted is the electrical transport converted to 
thermal units using the Wiedemann Franz law (Lq / p) . 
The two data sets lie on top of one another indicating 
that the WF law is satisfied as it should be in a regular 
metal. 

At intermediate fields, the thermal conductivity rises 
up to meet the normal state in a striking fashion. Fig. 
[3] shows the T — > extrapolated values as a function 
of magnetic field, where an applied field of only 20 mT 
doubles kq/T. We note that n/T continues to evolve with 
field up to 700mT, and remains unchanged above. This 



allows a convenient definition of H C 2 = 700 mT from a 
bulk transport probe. This is slightly below the value of 
H c2 = 900 mT reported in the previous study [20] . 

In Fig. [3] we also show the temperature dependence of 
the conductivity at T = 0.46 K to compare to Yamashita 
data. The three curves are all quantitatively consistent 
having an initial steep increase, followed by a plateauing 
of the conductivity value before increasing again towards 
the normal state value. The features are most appar- 
ent in the extrapolated zero-temperature data where the 
contribution is entirely electronic. The finite tempera- 
ture data are contaminated with phonons and it is clear 
from Fig. [I] that were the field dependence considered 
at a slightly higher temperature (T > 0.5 K) , the curve 
would be qualitatively different since the initial change 
would be a decrease in the conductivity. 

In a multi-band theory, there is little difference be- 
tween thermal conductivity of a d-wave or nodal s± sym- 
metry superconductor. In the low field regime, the nodal 
quasiparticles couple to the superfluid flow around a vor- 
tex and the energy states are doppler-shifted giving rise 
to a ~ y/H behaviour [37]. As the magnetic field is in- 
creased, the fully-gapped band has its gap suppressed 
and produces an additional channel for conductivity that 
causes the secondary upturn as H C 2 is approached. The 
low field doppler-shift behaviour has been observed in 
cuprate superconductors [231 [35] , but since these are es- 
sentially single band materials, the higher-field physics is 
not relevant. 

In order to make a qualitative comparison with the 
theoretical magnetic field dependence of the conductivity 
reported in I27j , we plot the normalised zero-temperature 
extrapolated values of the thermal conductivity as a func- 
tion of normalised field in Fig. [3j There is certainly a 
high degree of similarity, although quantitatively, the ra- 
tio of the magnitude of the increase in conductivity ef- 
fect in each of the different regimes does not appear to 
be consistent. Again, tuning the various parameters may 
allow for a better, possibly quantitative, agreement to 
be found. Beyond this, the magnetic field data supports 
both d-wsve and nodal s± symmetries. 

CONCLUSION 

In conclusion we have studied the temperature and 
field dependence of thermal transport in samples of 
LaFePO to temperatures as low as 60 mK. We find unam- 
biguous evidence for low energy excitations, interpreted 
to arise from nodal quasiparticle excitations. A quantita- 
tive comparison to theoretical models of heat transport in 
the pnictide superconductors suggest our data is consis- 
tent with an s± nodal scenario, which better captures the 
linear in T dependence of the zero field electronic ther- 
mal conductivity K e /T, and the insensitivity of T c to the 
measured large normal state scattering rate. Although 
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we cannot strictly rule out a d-wave symmetry, we find 
no quantitative agreement with scattering models in ci- 
ther the Born or unitary limit. We suggest that a study 
of universality in LaFePO could help further distinguish 
between the two scenarios. 
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